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ABSTRACT 


Models for the origin of the slow solar wind must account for two seemingly- 
contradictory observations: The slow wind has the composition of the closed-field 
corona, implying that it originates from the continuous opening and closing of 
flux at the boundary between open and dosed field. On the other hand, the slow 
wind has large angular width, up to ~ 60°, suggesting that its source extends far 
from the open-closed boundary. We propose a model that can explain both obser- 
vations. The key idea is that the source of the slow wind at the Sun is a network 
of narrow (possibly singular) open-field corridors that map to a web of separatri- 
ces and quasi-separatrix layers in the heliosphere. We compute analytically the 
topology of an open-field corridor and show that it produces a quasi-separatrix 
layer in the heliosphere that extends to angles far from the heliospheric current 
sheet. We then use an MHD code and MDI/SOHO observations of the photo- 
spheric magnetic field to calculate numerically, with high spatial resolution, the 
quasi-steady solar wind and magnetic field for a time period preceding the Au- 
gust 1, 2008 total solar eclipse. Our numerical results imply that, at least for this 
time period, a web of separatrices (which we term an S-web) forms with sufficient 
density and extent in the heliosphere to account for the observed properties of 
the slow wind. We discuss the implications of our S-web model for the structure 
and dynamics of the corona and heliosphere, and propose further tests of the 
model. 


Subject headings: Sun: magnetic field — Sun: corona 


Sun: solar wind 



1. Introduction 


Decades of in situ measurements of the heliosphere have firmly established that the Sun's 
wind consists of two distinct types: “fast” and “slow”. In terms of its origins at the Sun, the 
best understood is the fast wind, which typically exhibits speeds in excess of 600 km/s at 1 
AU and beyond (e.g., McComas et al. 2008). The fast wind is measured to be approximately 
steady, except for some Alfvenic turbulence. This wind is known to originate from coronal 
holes, regions that appear dark in XUV and X-ray images, due to a plasma density that is 
substantially lower (< 50%) than in surrounding coronal regions (Zirker 1977). As implied 
by eclipse and coronagraph images, the magnetic field in coronal holes is open appearing 
mainly radial and stretching out without end — whereas the field in the surrounding regions 
is closed, looping back down to the photosphere. Hence, the fast wind corresponds to the 
steady wind predicted by Parker in his classic work (Parker 1958, 1963). 

The slow wind, however, is much less understood. In particular, its origin at the Sun 
has long been one of the major unsolved problems in solar /heliospheric physics. This wind 
has a number of observed features that distinguish it physically from the fast wind. First, 
its speeds are typically slower, < 500 km/s. More important, it appears to be inherently 
non-steady. It exhibits strong and continuous variability in both plasma (for example, speed 
and composition) and magnetic field properties; variability that cannot be described as 
simply Alfvenic disturbances superimposed on a steady background. Finally, its location in 
the heliosphere is distinct; it is generally found surrounding the heliospheric current sheet 
(HCS) (e.g., Burlaga et al. 2002). A key point is that the HCS is always embedded inside 
slow wind, never fast. From the presently available spacecraft observations, it is not possible 
to rule out the possibility that slow wind also occurs in locations unconnected to the HCS, 
in other words, that there are pockets of slow wind with no embedded HCS and surrounded 
completely by fast wind. However, the present data are certainly consistent with the picture 
that all slow wind connects topologically to the HCS. 

Another key feature of the slow wind is its latitudinal extent, which typically ranges 
from 40° 60° near solar minimum, a time when it is easiest to distinguish the sources of fast 
and slow wind. Within this broad region of slow wind the actual HCS, across which the 
magnetic field changes direction, is very narrow. As for any current sheet, one can identify 
in the heliospheric data a scale over which the field becomes small and the plasma beta 
becomes large. This region is termed the plasma sheet and is usually of the order of a few 
degrees in angular width (e.g., Winterhalter et al. 1994; Bavassano et al. 1997; Wang et al. 
2000; Crooker et al. 2004). It is important to note that the HCS is often not symmetrically 
located within the broad band of slow wind, but is often found nearer to one edge of the 
slow wind region (Burlaga et al. 2002). It is also important to note that the field almost 
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never vanishes at the HCS, as would he expected for a true steady-state. This observation 
implies that, at least, the wind near the HCS must be continuously dynamic. 

The final and most critical feature of the slow wind that distinguishes it from the fast 
is the plasma composition (Geiss et al. 1995; von Steiger et, al. 1995). It is well-known that 
in the closed field corona, the ratio of the abundances of elements with low first ionization 
potential (FIP), such as Mg and Fe, to those with high FIP, such as C and Ne, is a factor 
4 or so higher than in the photosphere (e.g., Meyer 1985; Feldman & Widing 2003). This 
so-called FTP effect is not seen in the fast wind, which has abundances similar to those of the 
photosphere; but, it is present in the slow wind, which has abundances similar to that of the 
closed corona (Zurbuchen & von Steiger 2006; Zurbuchen 2007). Along with the difference in 
elemental abundances, the slow and fast wind also exhibit clear differences in their ion charge 
state abundances, for example, the ratio of 0 7 /0 6 . This ratio can be used to determine the 
“freeze- in” temperature of the ions at the source of the wind. The data show that the slow 
wind has a higher freeze-in temperature (> 1.5 x 10 6 K) than the fast wind (< 1.2 x 10 6 K) 
(von Steiger et al. 1997, 2001; Zurbuchen et al. 1999, 2002) 

The elemental abundances track very well the ionic abundances, indicating that there is 
a consistent compositional distinction between the two winds. Furthermore, the two winds 
have markedly different temporal variability in elemental and ionic composition. The fast 
exhibits an approximately constant composition: whereas the slow exhibits large and contin- 
uous variability, so that its elemental composition varies from coronal to near photospheric. 
The composition results suggest that the fast wind has a unique origin, presumably in coronal 
holes, but that the slow wind originates from a mixture of sources. 

In fact, Zurbuchen and coworkers have argued that the compositional differences, rather 
than the speed, are what truly distinguish the two winds, because it is possible to find solar 
wind whose composition and constancy match that of the “fast wind,” but that has relatively 
slow speed, < 500 km/s (Zhao et al. 2009). Note also that, as determined by the composition 
measurements (Zurbuchen et al. 1999), the boundary between the slow and fast wind in the 
heliosphere is sharp, of order a few degrees in angular extent, much smaller than the angular 
width of the slow wind region, but comparable to that of the plasma sheet. We conclude, 
therefore, that the fast and slow winds are far more appropriately described as the steady 
and unsteady winds. 

Since the differences in plasma composition of the two winds must be due to differences 
in their origins at the Sun, the composition data place severe constraints on the possible 
sources of the slow wind. In particular, the data imply that the slow wind originates in the 
dynamic opening of closed magnetic flux, which releases closed-corona plasma into the wind. 
Such a process would also naturally explain the difference in variability between the fast and 



It should be emphasized, however, that this constraint on the slow wind’s origin is 
not universally accepted. Several authors have argued that the slow wind originates from 
open-field coronal holes, just like the fast wind, but from the edges of the holes, where 
the field expands super-radially as it extends from the photosphere out to the heliosphere 
(e.g., Kovalenko 1981; Wang & Sheeley 1991: Cranmer k van Ballegooijen 2005; Cranmer 
et al. 2007). The hypothesis is that a large expansion factor can both slow down the wind 
by affecting the location of wave energy deposition in coronal flux tubes, and change the 
plasma composition by the FIP mechanism proposed by Laming (2004). The difficulty 
faced by the expansion factor models is that observations indicate that wind speed is not 
tightly correlated with composition. The wind from small equatorial coronal holes with a 
large expansion factor is indeed slow, with speeds < 500 km/s, in good agreement with the 
predictions of the expansion factor models. But this wind has photospheric FIP ratios, so 
it is still considered to be “fast wind” (Zhao et al. 2009). Furthermore, this not-so-fast 
wind has the temporal quasi-steadiness of the fast wind, rather than the quasi-chaotic time 
variation of the slow wind. 

We conclude, therefore, that the most likely source for the true slow wind, that with 
F IP-enhanced coronal composition, is the closed-field corona. In this case, the process that 
releases the coronal plasma to the wind must be either the opening of closed flux or inter- 
change reconnection between open and closed magnetic fields lines. This latter process is 
the underlying mechanism invoked by Fisk and co-workers (Fisk et al. 1998; Fisk 2003; Fisk 
k Zhao 2009) in their model for the heliospheric field. These authors argue that open flux 
can diffuse freely throughout the solar surface, even deep inside the helmet streamer region. 
If so, then the interchange reconnection between open and closed magnetic field lines would 
naturally account for both the composition and geometrical properties of the slow wind. 
The difficulty with this model is that it has not been demonstrated conclusively that this 
can actually occur. In fact, detailed MHD simulations indicate that it is difficult to bring 
open fields into closed-field regions without having them close down (Edmondson et al. 2010: 
Linker et al. 2010). These simulation results are in agreement with Antiochos et al. (2007), 
who argued that for the low-beta corona basic MHD force balance forbids the presence of 
open flux deep inside the closed helmet streamer region. 

Within the context of MHD models, the most likely location for the release of closed- 
field plasma is from the tops of helmet streamers (the Y-point at the bottom of the HCS), 
where the balance between gas pressure and magnetic pressure is most sensitive to perturba- 
tions. A number of authors have argued that streamer tops are unstable and should undergo 
continual opening and closing as a result of thermal instability (Suess et al. 199G; Endeve et 



al. 2004; Rappazzo et al. 2005). Even if streamer tops are stable, it seems inevitable that 
the constant emergence and disappearance of photospheric flux and the constant motions of 
the photospheric would force them to be continuously evolving. Furthermore, eoronagraph 
observations often show the ejection of “blobs” from the tops of streamers and into the HCS 
(Sheeley 1997). 

Although this streamer top model seems promising in that it naturally explains both 
the composition and variability, it has difficulty in accounting for the large angular widths 
of the slow wind. One would expect the instabilities to be confined to the high-plasma beta 
region about the current sheet. In fact, the plasma emanating from the streamer tops, the 
so-called stalks, is observed to be only ~ 3° 6 C wide, which agrees well with the plasma sheet 
width in the heliosphere (Bavassano et al. 1997; Wang et al. 2000). Even if the plasma sheet 
width were to be widened by the Kelvin-Helmholtz instability (e.g., Einaudi et al. 1999), 
there almost certainly would not be enough mass flux from the narrow region at the streamer 
tops to account for the slow wind. The streamer-top models can account for a thin band of 
slow wind around the HCS, but it seems unlikely that this is the origin of the bulk of the 
slow wind, which can extend as far as 30° in latitude from the HCS. 

In order to be compatible with the in situ data, we require some process that releases 
closed-field plasma onto open field lines that, in the heliosphere, can be far from the HCS. 
This requirement seems impossible to satisfy, because the plasma release must occur at the 
boundary between the open and closed field in the corona, which maps directly to the HCS. 
We describe below, however, a magnetic topology that resolves this slow wind paradox: the 
flux associated with an open-field corridor can be simultaneously near to and far from the 
open-closed boundary! 


2. The Topology of an Open-Field Corridor 

Figure 1 illustrates the magnetic connectivity from the photosphere to the heliosphere 

that results from an open-field corridor. The dark yellow inner sphere in the figure represents 

the photosphere, while the light yellow, semi-transparent one represents an arbitrary radial 

surface in the open-field heliosphere, say at 5/7 . The green line on the photosphere marks 

the boundary between open (gray) and closed (yellow) field regions, which is mapped by the 
magnetic field (red lines) to the HCS (thick green line) at the 5 R e surface. The green line 
at the HCS is also the polarity inversion line at this surface. Note that the four points, a, 
b, e, and d, which are meant to represent the end-points of the corridor at the Sun, map 
sequentially to the corresponding points a 2 * * * * 7 , 1/. c'. and d' along the HCS. 
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The open field pattern at the photosphere of Fig. 1 consists of a large polar coronal 
hole and. as is often seen, a smaller low-latitude hole. In recent work, we argued that if the 
two holes are in the same photospheric polarity region, then by our uniqueness conjecture 
the holes must be connected by an open field corridor, as illustrated above (Antiochos et al. 
2007). It is evident from the figure that the flux in the corridor maps on the heliospheric 
surface to a thin arc (light, gray band), bounded at both ends by the HCS. The flux between 
the arc and the HCS maps to the low-latitude extension while the flux outside the arc maps 
to the main part of the polar coronal hole. The corridor and its associated arc are the 
footprints of two quasi-separatrix layers (QSLs, e.g., Priest k Demoulin 1995; Demoulin et 
al. 1996) that combine into a hyperbolic flux tube, as has been described in detail by Titov 
et al. (2002, 2008) for the case of closed magnetic configurations. In contrast, the HCS is a 
true separatrix. 

The key point for understanding the origin of the slow wind is that, just like the HCS, the 
QSL arc in the heliosphere is also a source region for slow wind. If the open-field corridor 
at the Sun is sufficiently narrow, then the continual evolution of the photosphere, driven 
by the ever-present supergranular flow and flux emergence/submergence in particular, will 
continually change the exact location of this corridor. But, by the uniqueness conjecture 
(Antiochos et al. 2007), the corridor is a topologically robust feature, similar to a null-point, 
and must be present on the photosphere as long as the low-latitude coronal hole extension 
is present. Its location and shape, however, will vary in response to local photospheric 
changes. These variations require field line opening/closing and interchange reconnection, 
thereby releasing closed-field plasma all along the QSL arc in the heliosphere. Therefore, if 
the QSL arc extends to high latitudes, this will naturally produce slow wind with an extent 
far from the HCS. 

To determine whether the QSL resulting from an open field corridor does, indeed, reach 
high heliospheric latitudes, we have calculated an example of a field such as that of Fig. 1 
using the source surface model (Altschuler k Newkirk 1969; Schatten et al. 1969; Hoeksema 
1991). The field is most easily determined from the image-dipole formula derived by Anti- 
ochos et al. (2007). For a dipole with moment d located at a point r d inside the Sun, and 
a source surface at radius R s , the magnetic field B is determined from the potential via 
B = — V<f>. where <f> is given by: 

... d - (r - r d ) Rsr'jd • (R%r - r 2 r d ) 

|r-r f/ | 3 \rjr-K 2 s r d \ 3 ' ' 

This field satisfies the source-surface boundary condition that B $ = I3 # = 0 at r = R s , 
since d> = 0 there. The advantage of this formulation is that most active regions can be 
approximated by a collection of dipoles, and one can build up a field of arbitrary complexity 
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by simply adding a series of dipoles of the form of Eq. (1). Each dipole is specified in terms 
of its position in spherical coordinates r d = r d r(O d , <p d ), where r d , 6 d , and <f> d specify the 
location of the dipole, and the spherical components of its dipole moment, d = (d r . cl#, d#). 

Figure 2 shows the field computed from Eq. (1) for the case of two dipoles: a sun-centered 
global dipole with a dipole moment of unit magnitude directed along the north polar axis, 
and an equatorial “active region” dipole at r d = O.8i? 0 r (90°. 0°) with a northward-pointing 
dipole moment d = (0, -0.2,0). The source surface radius is chosen as R s = 4/? . . though 
the exact value is not critical for our argument. Note that for convenience in viewing the 
magnetic field, we have selected the dipole parameters so that the system has symmetry 
across both the equatorial ( 6 = 90°) and meridional {<i> = 0) planes. Also, for ease of 
viewing, we show in the Fig. 2 only the front hemisphere defined by the angular region 
(15° < 9 < 90°) and (-90° < </> < 90°). 

The solar surface, the photosphere, corresponds to the gray grid in Fig. 2. The colored 
contours on this surface correspond to contours of radial flux, indicating the presence of the 
active region dipole at the equator. We selected the parameters for the active region dipole 
so that its structure would be easily resolved. It is evident from Fig. 2 that the region is large 
compared to real active regions, which are generally only a few degrees in angular extent. 
On the other hand, the maximum field strength at the dipole center is only ~ 20 times that 
of the polar region, which is much less than the corresponding ratio for solar active regions, 
so the flux ratio between the active region and global background field is approximately 
correct. This ratio is the important parameter to obtain a. coronal hole extension. 

The thick black line along the equator is the B r = 0 contour, i.e., the polarity inversion 
line. The thick black line above the solar surface is the polarity inversion line at the source 
surface, i.e., the bottom of the HCS. Red field lines are traced at equal intervals along the 
HCS down to the solar surface. These define the boundary between open and closed field 
lines. As expected, the effect of the equatorial dipole is to pull the open-closed boundary 
down to lower latitudes; in other words, to create a low-latitude extension of the coronal 
hole, which can be seen as the gray shaded region in the Figure. Far from the dipole, the 
coronal hole boundary is at a latitude of ~ 54°, whereas at the meridional symmetry plane 
the boundary drops down to ~ 26°. 

For the large spatial scale of our active region dipole, the extension of the coronal hole 
down to low latitudes is gradual rather than in the form of a distinct “elephant trunk” , but 
the basic effect is clearly present. There is no open-field corridor in Fig. 2, but let us now add 
another dipole to the system, displaced 20° in both latitude and longitude from the equatorial 
one and a factor of five times weaker. This dipole is located at r d = 0.8f? e f(70°, 20°) 
with a primarily southward-pointing dipole moment d = (0,0.05,0). In order to maintain 



the equatorial and meridional symmetry, as mentioned earlier, we actually add 4 dipoles 
symmetrically located about the equatorial and meridional planes. 

The resulting field is shown in Figure 3. The effect of the additional dipoles is to add 
high-latitude polarity inversion lines to the system. These “squeeze” the open-flux extension 
of Fig. 2 to form a narrow corridor and a low-latitude coronal hole. As in Fig. 2, red field 
lines are traced from equidistant foot points along the HCS down to the solar surface. The 
red footpoints at the photosphere appear to traverse the boundary of the low-latitude hole 
and then jump abruptly to the polar hole boundary, which implies that the mapping defined 
by the field develops extreme gradients in the region connecting the two holes. To clarify 
this point, we have traced two sets of field lines, colored in blue, from footpoints that are 
closely located at the HCS. The corresponding solar footpoints are much more widely spaced, 
running along the corridor. The resulting structure, Fig. 3, looks very similar to the mapping 
drawn in Fig. 1, in that the closely spaced pairs of points abb' and c',d' at the HCS map to 
far-separated points a.b and c.d at the solar surface. Note also that although the footpoints 
of the two sets of blue lines approach each other very closely at the photosphere, they are 
far separated at the HCS, by a distance of order i? 0 . This result indicates that even though 
the low-latitude coronal hole has small area, it contains a substantial magnetic flux. As is 
evident from the colored contours in Fig. 3, the photospheric field strength in the low-latitude 
hole is large due to the presence of the active region dipole. 

The analytic model underlying Fig. 3 has similar topology to the case shown schemat- 
ically in Fig. 1. The low-latitude coronal hole extension in Fig. 3 is connected to the main 
polar hole by a corridor that becomes very narrow. Furthermore, this type of topology is not 
difficult to obtain. It is often observed in quasi-steady MHD solutions for observed photo- 
spheric fields, as will be shown below. A similar corridor was found for Carrington rotation 
1922 (Antiochos et al. 2007). 

The question now is whether the open flux in the corridor connects to large latitudes 
in the heliosphere. To answer this question, we trace field lines from a, set of photospheric 
footpoints lying on a latitudinal line segment spanning the narrowest width of the corridor, 
which is only of order 5.000 km at the photosphere. Fig. 4b shows the footpoints and the 
field lines (green) near the photosphere and Fig. 4a shows where they map to on the source 
surface. We note that the corridor maps to high latitudes. In fact, for this analytic case, the 
corridor mapping defines a QSL arc that reaches latitudes > 45°, greater than that of the 
observed slow wind. 

This result, that the corridor maps to heliospheric latitudes far above the HCS, is 
robust in that it is not sensitive to the exact position of the secondary dipole. The position 
and geometry of the corridor, on the other hand, is very sensitive to the photospheric flux 



distribution. For example, its width would change and its location would change substantially 
if the secondary dipole were moved in longitude. Based on flux conservation arguments, and 
the fact that the heliospheric magnetic field is almost uniform in latitude, we can argue that 
the angular extent of the QSL arc, however, would be expected to depend primarily on the 
ratio of the flux in the low-latitude coronal hole extension to that in the polar hole. For 
example, in the extreme case that the fluxes were equal, the corridor mapping would be 
expected to reach the heliospheric pole (90° from the HCS!), irrespective of the geometry of 
the corridor or of the coronal holes. 


3. The S-Web Model 

If the width of the corridor at the photosphere is small compared to the scale of typical 
motions there, such as the supergranular flow, we expect that the whole corridor will contin- 
uously disrupt and reform at the photosphere and, consequently, closed- field plasma will be 
released by reconnection all along the QSL arc in the heliosphere. Therefore, the topology 
of Fig. 2 can resolve the slow wind paradox. The overriding question, however, is whether 
there are enough such corridors and corresponding QSL arcs in the heliosphere to account 
for the slow wind that is observed. The' flux distribution of Fig. 2 produces only one such 
arc, which would likely not be sufficient to reproduce the observed slow' wind. 

In order to answer this question, we calculated the quasi-steady model for an observed 
photospheric flux distribution. Figure 5a shows the photospheric radial field as derived from 
MDI observations on SOHO (Scherrer et al. 1995) for a time period preceding the August 
1, 2008 total solar eclipse. This calculation was used to predict the structure of the corona 
prior to the eclipse, using magnetic field data measured during the period June 25- July 21, 
2008. The prediction compares very favorably with images of the corona taken during the 
eclipse in Mongolia (Rusin et al. 2010). Note that the high resolution of the calculation 
captures the details of many small-scale bipoles in the photospheric magnetic field (Harvey 
1985). This has been incorporated into the idea of the “magnetic carpet” (Schrijver et al. 
1997). We also show the polarity inversion line B r = 0 slightly above the photosphere, at 
r = 1.05 i?© to delineate the magnetic polarity of the large-scale structures. (The polarity 
inversion line in the photosphere itself shows an enormous complexity that overshadows its 
usefulness to discern the large-scale magnetic polarity.) 

The quasi-steadv model was calculated by using the 3D MHD code MAS. The MAS 
code and its implementation are described in detail by Mikic & Linker (1994), Mikic et al. 
(1999), Linker et al. (1999), and Lionello et al. (2009). MAS solves the time-dependent 
MHD equations, including a realistic energy equation with optically thin radiation and ther- 



mal conduction parallel to the magnetic field. Given the magnetic field at the photosphere 
and an assumption for the coronal heating source, the MHD equations are advanced until 
the magnetic field settles down close to steady state. MHD models are generally considered 
to be the most sophisticated implementation of Parker’s solar wind theory because they 
incorporate all the essential physics, including the balance between gas pressure and Lorentz 
force. An important assumption is the form of coronal heating, which is prescribed empiri- 
cally at the present time since the coronal heating process is still unknown. The parameters 
of the empirical heating model are constrained by observations of coronal emission in EUV 
and X-rays (e.g., Lionello et al. 2009), as well as by solar wind measurements. 

In order to capture as much of the photospheric magnetic structure as possible, we ran 
the MAS code with unprecedented resolution. Our calculation used more than 16 million 
mesh cells and was run on over 4000 processors of NSF’s Ranger supercomputer at the Texas 
Advanced Computing Center, making it possible to include much of the small-scale structure 
of the photospheric field in both the quiet sun and in coronal holes, as shown in Fig. 5a. 
These calculations are unique in the degree to which they capture the small-scale structure 
of the measured magnetic field. 

Figure 5b shows the distribution of open and closed magnetic field regions at the solar 
surface as determined by the model. It is evident that there are many low-latitude coronal 
hole extensions, similar to that in Fig. 3, but with much more structure. Several of these 
extensions appear to be disconnected from the main polar holes, but this is partly due to the 
limited resolution of the figure. A few of these coronal hole extensions are indeed connected 
by very thin corridors in the photosphere, though many are only linked to the polar coronal 
holes in a singular manner, as described in detail by Titov et al. (2010), and as discussed 
further below. 

Although the open field pattern in Fig. 5b is very complex, the important issue is the 
degree of complexity of the open field that maps into the heliosphere and, in particular, the 
structure of the separatrices and QSLs there. We determined the open field mapping in great 
detail by tracing tens of millions of magnetic field lines. The topology of this mapping, as 
evidenced by structures such as separatrices and QSLs, is most easily seen by analyzing the 
squashing factor Q (Titov et al. 2002; Titov 2007). Q is a measure of the distortion in the 
magnetic field mapping, and is directly related to the gradients in the connectivity. QSLs 
are regions of very large Q: we generally define them as any region with Q > 10 3 . True 
separatrices such as the HCS have infinite Q, because the mapping is singular there, but 
when computed numerically they appear as surfaces with very large (unresolved) values of 
Q. The gray arc* at r — 5i? e hi Fig. 1 is a QSL in the open field, and consequently would be a 
region of high Q. The green HCS would also be a region of high (infinite) Q. As will be seen 



11 - 

below, a high-resolution analysis of the Q properties of our MHD simulation is extremely 
informative. 

Figure 6a shows Q in a meridional plane at a central Carrington longitude of 23.33° at 
the time of the eclipse at 10:21UT, while Figure 6b shows magnetic field lines traced from the 
vicinity of the solar limbs at the same time. We see that Q outlines the boundary between 
open and closed field, which is a true separatrix surface, but it is apparent that there is much 
more detailed structure in both the closed and open field regions. The complex structure 
of Q in the closed-field region is expected; it simply reflects the fact that the photospheric 
field consists of many small bipoles; but, there is also substantial structure in the open field 
near the open-closed boundary. Note the presence of a “pseudostreamer" on the NE limb, a 
feature that has been discussed by Wang et al. (2007). The relationship of pseudostreamers 
to open hole corridors and the S-web is discussed in detail in Titov et al. (2010) 

Figure 7a shows Q in the spherical surface at r = 10/?® using a logarithmic scale. This is 
the structure that is expected to map into the inner heliosphere (appropriately wrapped into 
a spiral magnetic field by solar rotation), since the magnetic field has reached its asymptotic 
structure by this radius. The thick black line is the heliospheric current sheet (at which B r 
reverses sign). Figure 7b shows the magnitude of B r at the same radial surface r = 10/?®. 
Note that the choice of 10/?® is not crucial. Any surface in the heliosphere (where the field 
is all open) yields similar results. 

It is important to emphasize that the apparent structure in Q expresses only the con- 
nectivity of the open field, not its actual magnitude. In spite of the enormous magnetic 
complexity at the solar surface, the radial field distribution in the heliosphere is completely 
unremarkable, Fig. 7b. There is a single polarity inversion line denoting a single HCS, as 
is generally observed near solar minimum, arid this HCS runs more or less equatorial. The 
radial field is essentially uniform away from the HCS, as would be expected from simple pres- 
sure balance. (Careful examination of the plot of B r shows that there is a faint semblance 
of the structure that can be seen in Q, but it is only a small perturbation.) 

On the other hand, the Q map at this surface is remarkable, indeed, Fig. 7a. We see 
that surrounding the HCS is a broad web of separatrices and QSLs of enormous complexity. 
There are at least four striking features of this S-web. First, it has an angular extent in 
latitude of approximately 40°. sufficient to account for the observed extent of the slow wind. 
Note also that the angular extent does vary with longitude, but only by a factor of two or 
so. Second, the HCS is not necessarily in the center of the S-web. but is sometimes near 
its edge. This can explain the frequent observation that the HCS is usually not centrally 
located within slow wind streams (e.g., Burlaga et al. 2002). Third, the boundary between 
the S-web layer and the featureless polar hole region is sharp; it is narrow compared to the 
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width of the S-web. This can explain the observation that the transition from slow to fast 
wind as measured by the composition data is narrow compared to the slow wind region itself 
(Zurbuchen et al. 1999). 

In order to explore the details of how coronal hole extensions connect to the polar holes, 
we calculated coronal hole areas at different heights in the corona. Figure 8 shows the location 
of a region near longitude 75° and latitude 15°N in which we explored the connection between 
the low-latitude coronal hole extensions (of negative polarity, shown in blue) in detail. It is 
evident that the coronal hole extensions in this region appear disconnected from the north 
polar hole in the photosphere, but connect with it low in the corona (at heights approximately 
between O.Oli?© and (1.02/7. above the photosphere). Figure 9 shows explicitly how these 
coronal holes connect in the low corona. The three-dimensional shape of the coronal hole 
boundary is shown as a green semi-transparent surface in the low corona in the region detailed 
in Figure 8. This is the boundary between open and closed field regions. The regions marked 
by A, B. and C show examples in which the extensions of coronal holes are not connected 
in the photosphere, at least by any measurable open-field corridor, but appear to connect 
above the photosphere in the low corona. These regions are also indicated in Figure 8 for 
ease of cross-reference. Despite the fact that these coronal holes are “disconnected” in the 
photosphere, they always remain topologically linked in a singular manner with the polar 
coronal hole, as discussed by Titov et al. (2010). 

Finally, note that the connections of the high-Q lines between the neighborhood of the 
HCS and the photosphere and low corona that were postulated by the uniqueness conjecture 
(Antiochos et al. 2007) are largely present, even though the insight from these new high- 
resolution MHD simulations has led us to generalize the uniqueness conjecture. We have 
found that, in general, coronal hole extensions are sometimes connected to the polar holes in 
the photosphere via narrow corridors, as originally postulated (Antiochos et al. 2007), but in 
other instances they are disconnected in the photosphere, but remain topologically linked to 
the polar holes (Titov et al. 2010). In either case, these connections are responsible for the 
formation of the S-web. It should be emphasized that discovery of the extremely complex 
structure of these connections would not have been possible without our high-resolution 
models and analysis tools. 


4. Discussion 


The major conclusion from our results is that the underlying premise of the streamer 
top model is correct. The slow wind is expected to originate from the release of closed- 
field plasma due to the dynamic rearrangement of the open-closed field boundary. The 



key new addition of our S-web model to this picture is that the inherent complexity of the 
photospheric field leads to a network of narrowly connected and disconnected coronal holes 
that nevertheless always remain linked. This produces a separatrix web in the heliosphere 
that extends the release of slow wind to regions that significantly depart from the HCS. 
Hence, our model accounts for both the observed composition and the extent of the slow 
wind. 

One immediate prediction from the model is that the angular width of the slow wind 
is determined primarily by the complexity of the flux distribution in the photosphere. This 
complexity produces a very convoluted polarity inversion line in the low corona and an intri- 
cate coronal hole pattern (Figure 5). Our ability to identify the S-web and its manifestations 
rests on high-resolution calculations that are beginning to capture the multitude of small 
dipoles in the photospheric magnetic field. If the solar field were a pure dipole, producing 
an inversion line that runs straight along the equator, then only the polar coronal holes 
would be present and there would be no separatrix web in the heliosphere. For this “basal” 
(though idealized) slow wind case, if we assume that the dynamic broadening of the open- 
closed boundary at the Sun is of order the scale of a supergranule, ~ 30,000 km, the angular 
extent of the wind would be only of order 3° 5°, and would be centered about the HCS. Of 
course, the solar field is never a simple dipole. 

At the present time we do not know if the complexity seen in Figures 5 7 is typical, 
or whether it is particular to this late declining phase of Cycle 23. It should be noted that 
the present minimum appears to be somewhat different than the previous few minima. In 
particular, the polar field strength is significantly weaker (e.g., Luhmarm et al. 2009). 

The S-web model predicts that for time periods during which extensions of coronal holes 
away from the main polar holes are leas extensive than in Cycle 23, the angular extent of 
the slow wind region would be smaller. In fact, there is evidence from radio scintillation 
data (Manoharan, P. K. 2009, private communication) and recent Ulysses solar wind mea- 
surements that the Cycle 23 minimum has a substantially broader slow wind region than 
that of the previous cycle. Indeed, during the previous minimum (circa 1996), equatorial 
coronal hole extensions were less common than during the recent solar minimum. Further 
high-resolution numerical calculations will be needed to address this result. 

Another prediction of the model is that the slow wind region is actually a mixture of 
winds. It is evident from Fig. 7 that the separatrix web is not space- filling. There are regions 
within the broad S-web band where the wind emanates from the low-latitude coronal hole 
extensions. These regions are likely to have large expansion factor, so that the wind will be 
slow compared to the fast wind from the polar regions, but its composition will be different 
than that of closed-field plasma. Our model, therefore, naturally explains the observed 



variability of the slow wind composition. 

A key aspect of the S-web model that has yet to be calculated is the dynamic release 
of closed-field plasma. Although our quasi-steady calculations allow us to investigate the 
topology of the field, and to identify the structure of the separatrix web in the heliosphere, 
they do not actually produce a slow wind with closed-field composition yet. For this we 
need fully dynamic simulations that, include the driving due to photospheric motions (e.g., 
resulting from differential rotation) and flux emergence. Such simulations are now being- 
performed in 3D (e.g., Edmondson et al. 2009, 2010; Linker et al. 2010) for simplified pho- 
tospheric flux distributions and driving flows. These simulations do verify the basic idea of 
the S-web model that open-field corridors will form and evolve in response to photospheric 
motions (Edmondson et al. 2009). Higher resolution simulations will be needed, however, 
to test the model in detail. On the other hand, it seems unlikely that dynamic calculations 
with the degree of structure present in Fig. 7 will lie feasible in the near future. It may well 
be that a definitive treatment of the slow wind will require the development of a statistical 
theory of the dynamics of the S-web model. 
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Fig. 1. — Magnetic field topology of an open field region consisting of a large polar coronal 
hole and a smaller low-latitude hole connected by an open-field corridor. The inner surface 
is the photosphere, with the dark gray and bright yellow regions corresponding to open and 
closed field respectively. The outer transparent surface is a radial surface in the heliosphere. 
The dark green line is the polarity inversion line and the light gray arc indicates where the 
open-field corridor maps to on this outer surface. 




Fig. 2. — (Top) Open-closed magnetic field topology for a photospheric flux distribution 
due to a global dipole and an equatorial dipole. The gray shaded region indicates the polar 
coronal hole (the open flux region). The contours on the inner surface indicate radial field 
magnitude at the photosphere. The black lines correspond to the polarity inversion line 
at the photosphere and source surface. The red lines are magnetic field lines. (Bottom) 
Close-up near the solar surface of the magnetic field above. 
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Fig. 5, (a) Distribution of the radial component of the magnetic field in the photosphere 

that was used in the MHD simulation to predict the structure of the corona for the August 1, 
2008 eclipse, as deduced from MDI measurements, (b) The open and closed field regions in 
the photosphere as determined from the MHD solution. The polarity inversion line (B r = 0) 
at a height r = 1.05/fy is superimposed on these images to aid in identifying the polarity of 
the large-scale magnetic flux. 
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August 1, 2008 Total Solar Eclipse 
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Fig. 6. (a) Plot of the squashing factor Q on a logarithmic scale in a meridional plane at 

the time of the eclipse on August 1, 2008 at 10:21UT. In this view, solar north is vertically 
up and the £ 0 angle is zero. [At the time of the eclipse B 0 = 5.8°, so this view is slightly 
different than what would have been observed.] The Sun’s surface is colored by the value of 
B r with the same scaling as that in Fig. 5. (b) Magnetic field lines traced from the vicinity 
of the limbs at the same time, showing the structure of the open and closed field regions. 





24 


August 1, 2008 Total Solar Eclipse 



Longitude 




90°N 


Radial Magnetic Field at r = 10R o 


60° 120” 180° 240“ 300” 360° 

Longitude 


-90°S 

0 ” 


Fig, 7. (a) Plot of the squashing factor Q in the spherical surface r = 10i?, 3 on a logarithmic 

scale versus longitude and latitude, (b) Plot of B r in the same spherical surface. The HCS 
(i.e.. the location of B r = 0) is superimposed on these images as a thick black line. The 
complex structure in Q in the vicinity of the HCS is produced by the S-web. 
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Fig. 8. The variation of coronal hole shape with height above the photosphere. The 
top panel shows coronal holes at the photosphere (the same information as in Fig. 5b) for 
reference. The black square shows a 100° x 100° region centered at longitude 75° and latitude 
15 X that was used to compute the detailed variation of coronal hole shape with radius in the 
lower panels. Note that the extended coronal holes connect to the polar holes low in the solar 
corona. The regions denoted by A. B, and C are cross-referenced with the corresponding 
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August 1, 2008 Total Solar Eclipse 


The Shape of the Coronal Hole Boundary 



Fig. 9. The three-dimensional shape of the coronal hole boundary (semi-transparent sur- 
face) in the region detailed in Figure 8, showing that some of the coronal hole extensions 
(blue areas on the surface of the sphere) connect with the north polar hole low in the corona. 
The top panel shows a view in which the surface is artificially stretched in radius by a factor 
of 3x to show details near the photosphere. The bottom left panel shows the same view 
without the radial stretching. The bottom right panel shows the region detailed in the con- 
text of the whole Sun. The regions denoted by A. B. and C are cross-referenced with the 
corresponding regions in Fig. 8. 





